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We investigate sound propagation in a moving fluid conflned in a randomly corrugated tube. For 
weak randomness and small fluid velocities the localization length ^ shows extreme sensitivity 
to the variation of In the opposite limit of large fluid velocities, (, acquires a constant value 
which is independent of the frequency of the incident sound wave, the degree of randomness and 
itself. Finally, we And that the standard deviation ai^T of the logarithm of transmittance ln(r) 
is a universal function of the ensemble average (In T) , which is not affected by the fluid velocity. 


Introduction - Wave transport in random media is the 
focus of many theoretical and experimental studies dur¬ 
ing the last sixty years [T]. One of the fascinating phe¬ 
nomena that has been predicted and subsequently ob¬ 
served in such media is the halt of wave propagation: it 
was found that due to multiple scattering and the conse¬ 
quent destructive interference between scattered waves, 
the total transmittance decays exponentially with the 
size of the system. This phenomenon, known as An¬ 
derson localization, has been originally predicted in the 
realm of condensed matter physics mm- More recently it 
has been studied and observed in optics gHI] , microwaves 
[SlIS] , acoustics m, as well as for matter waves in cold 
atoms systems mm- 

In the present paper we investigate Anderson localiza¬ 
tion in a new setting, namely, sound propagation in a 
moving fluid. For simplicity we assume that the fluid is 
confined to a one-channel waveguide with random cor¬ 
rugation. We consider subsonic flows where the fluid is 
inviscid and turbulent effects are irrelevant. We find that 
wave interference that lead to Anderson localization of 
sound are strongly affected by the velocity of the flow 
In a broad range of parameter (strength of disorder, 
wave frequency, size and cross-section of the scatterers) 
the localization length ^ is extremely sensitive to as 
long as is not too large. As increases, ^ saturates 
at a universal value, independent of the wave frequency. 
We also find that the variance of the logarithm of 
transmittance ln(T) is a universal function of its average 
value (InT), independent of the flow velocity 

Mathematical modeling - We consider sound propaga¬ 
tion in a tube with a moving fluid. Moreover we will focus 
our analysis on one- dimensional propagation which im¬ 
poses the constraint that the wavelength of the sound is 
much larger than the typical width of the tube. The tube 
consists of two parts: the left x < —L/2 and the right 
X > L/2 (semi-infinite) domains have a constant cross 
section Aq and constitute the “leads” from where the 
sound is emitted and detected respectively, while in the 
domain —L/2 < x < L/2 the tube cross-section A(x) is 
non-uniform (corrugated domain). We will assume that 


the fluid flows from left to right with a constant velocity 
at the leads. 

The one-dimensional equations for mass and momen¬ 
tum conservation read m- 

A(x) -f ^ [p{x, t)v{x, t)A{x)] =0; (1) 

dv{x,t) dv{x,t) 1 dp{x,t) 


where p{x, t) denotes the fluid density, v{x, t) is its veloc¬ 
ity, and p{x, t) is the pressure. In what follows we assume 
that the corrugated region consists of uniform segments, 
with cross-section A„ for the n-th segment. The seg¬ 
ments are separated by sharp transition regions in which 
the cross-section rapidly changes from An to An+i- The 
precise profile of various quantities in the transition re¬ 
gions are complicated and we eliminate those regions by 
imposing boundary conditions which are obtained by in¬ 
tegrating Eqs. Q across the transition regions between 
adjacent segments. This results in the continuity of the 
quantities: pvA and -f- w, where w is the enthalpy 
per unit mass. While integrating the second equation in 
Eq. ,0 we have used the relation {dw)s = (^)dp which 
is valid for isentropic flow (the entropy s per unit mass 
is constant). 

All the quantities in Eq. Q , for each segment n, have a 
stationary (time-independent) part upon 

which small oscillatory terms {p'{x,t), p'{x,t),v'{x,t)} 
are superimposed, due to the sound wave. Linearizing 
Eqs. ([^ with respect to the oscillatory terms yields, in 
the n—th segment 


Dp'n 

Dt 




dv^n 

dx 


= 0 ; 


Dv’n ^ PI 

Dt p(0) Qj. 


( 2 ) 


where = (^§1 + Vn ^is the so-called convective 

derivative and cq = \/ {dp/dp) a the adiabatic speed of 
sound in the reference frame of the fluid. In deriving 
Eqs. ([^ we have used the relation p/ = Cpp/. Moreover, 
to somewhat simplify the treatment, we have assumed 
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a nearly incompressible fluid so that cq is the same in 
each segment. Eliminating v'^ from Eqs. ([^ leads to the 
equation for the pressure wave in a moving fluid [14j 


^^p'n _ 2 9^Pn 
Dt^ ° dx^ ’ 


(3) 


whose solution can be written in terms of two counter- 
propagating waves: 


Pnix,t) 


jZif ik^x—iu)t 



x—iujt 


(4) 


where w = fc^(co -I- Vn^) = k^(co - Vn^) and kl{k\) in¬ 
dicate the wave vector of the waves propagating in (op¬ 
posite to) the direction of the flow. Correspondingly, the 
velocity variation v'^{x,t) as determined from Eq. (§ is 


At the boundary between two adjacent segments the 
values of v' and p' should be matched by the bound¬ 
ary conditions. The latter are obtained by lineariza¬ 
tion of the two aforementioned continuity conditions, 
namely for pvA and for 4 ^ 1 ^ -(-■{/;. This results in conti¬ 
nuity, accross the boundary, of the stationary quantities 
and -I- Wq, and of the oscillatory quanti¬ 
ties p''°\v'A) + {v^^^A)p', + -^pk The last two 

conditions can be conveniently re-written as 


j(0)^ 




1 - 



AnVn 


,(0) 


Pn+l 


■ V, 


(0) / 
n-l-l^n+1) 


( 6 ) 





At the left and right leads, where the cross-section Aq is 
constant, the pressure and velocity variations are given 
by Eqs. Q and ([^ with the sub-indexes n = L{R) for the 
left (right) lead, and . In the corrugated 

domain —L/2 < x < L/2, substitution of Eqs. Q and ([^ 
into Eqs. (|^ allows us to cast the boundary relations into 
a transfer matrix form which connects the forward and 
backward propagating sound wave amplitudes between 
the two subsequent domains, n and n -I-1. 

Finally, the pressure at a: = L/2 is related to that at 
X = —L/2 via the total transfer matrix Ad: 

/ pf pik^ L/2 \ ( pf pikA-L/2) \ 

y p^g-ik>’L/2 j ^ y pb^-ik\-L/2) j ■ 

The transmission and reflection amplitudes for left and 
right incident waves can be expressed in terms of the 
transfer matrix elements as tj, = ; ll = — 

~ relations have been ob¬ 

tained from Eqs. Q by imposing the appropriate scatter¬ 
ing conditions = 0 {P[ = 0) associated to left (right) 


incident waves respectively. Furthermore, straightfor¬ 
ward calculations show that |detAd| = 1 , which means 
that \tL\ = \tR\ = \t\. 

Transmittance T and reflectance R are defined as ra¬ 
tios of the corresponding energy fluxes (currents), I [15] . 
The latter are proportional to the square of the wave am¬ 
plitudes, for instance, the current for the forward propa¬ 
gating wave in the left lead is ^ However, the 

proportionality coefficients for the forward and backward 
propagating waves are different. Since transmittance in¬ 
volves a pair of waves (incident and transmitted) propa¬ 
gating in the same direction, one simply has Tr = 
for the transmittance from left to right and, similarly, 
Tr = for transmittance in the opposite direction 

(the two are the same due to the previously mentioned 
equality \tL\ = |tfl|). On the other hand, reflectance in¬ 
volves a pair of waves (incident and reflected) propagat¬ 
ing in the opposite directions which leads to the relations 
Rl = 7 • \ll\^ and Rr = (I/ 7 ) • \t'r\^- The coefficient 
/ _ (0) \ ^ 

7 = ( e°_)_ 0 (o) ) is due to the fact that the reflected waves 
travel with a different speed than the incident waves. 
Armed with the above knowledge we are now ready to 
investigate the effects of interference due to scattering 
from defects. 



FIG. 1: Pressure profile p'( 2 :) along a periodically modulated 
tube with corrugated areas characterized by Li — L 2 , Ai = 
O. 6 A 2 and a defect Ld = 0.051/2, Ad = A 2 in the middle of 
the tube (we consider L 2 = 1 ,A 2 = 1). A localized mode 
is created due to interference effects between the periodic 
modulation and the defect. The cross section of the leads 
is Ao = O. 6 A 2 and the velocity of the fluid is = O.Icq. 

One Defect - It is instructive to start our analysis with 
the simple example of one defect embedded in an other¬ 
wise uniform infinite tube of cross section Ag. The cross 
section of the defect is Ad and it occupies the interval 
—Ldl2 < X < Ld/2. First we note that, regardless of 
the ratio AdjA^, a single defect cannot support a bound 
state. This is because Eq. (§ in the leads, away from 
the defect, cannot have exponentially decaying solutions 
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for real w, i.e. the dispersion relation requires real k for 
real to. The situation here is different from that in quan¬ 
tum mechanics, where an attractive potential can pro¬ 
duce a bound state with a negative energy (imaginary 
k). However when the defect is placed in a periodically 
modulated tube consisting of segments with lengths Li, 
L 2 and cross-sections Ai, A 2 then a localized defect state 
can be formed (see Fig. [^. This is a result of interference 
between the periodic corrugation and the defect. 

Next we study the transmittance properties of a single 
defect. The boundary conditions Eqs. Q at a: = ±Ld/2 
allow us to evaluate the transfer matrix and conse¬ 
quently the transmittance T^. We obtain: 


that for small velocities < 0.4co, the localization 
length is sensitive to the variations of while for larger 
values of it originally oscillates and finally saturates 
to a universal value which is independent of w, 61 , 62 , and 
v^°'> itself. 

An understanding of the universal value of the rescaled 
localization length is achieved by employing a random 
phase approximation (RPA). We consider a tube consist¬ 
ing of {N — 1) scattering units to which we add one more 
scattering unit. The transmission amplitude of the 
combined system is 





( 10 ) 


where a = ^, (j)d = and 13d = v^f^/co- From 

Eq. ([^ we see that the resonance modes of the defect 
(corresponding to = 1) are achieved when ipd = 27rm 
(where m = 1, 2, • • •). The resonance condition can be re¬ 
written in a more transparent way as (fc;J -I- k^)Ld = 27rm 
which resembles the standard resonance condition, albeit 
now the two counter-propagating waves have different 
wave-vectors. Finally, it is important to realize that (pd 
depends not only on the incident frequency w but also 
on the velocity of the flow Pd which can lead to strong 
changes in transmittance. 

Disorder Case - We proceed with the analysis of the 
transport properties of the sound in a disordered tube. 
We therefore destroy the periodicity of the corrugated 
tube consisting of two different cross sections Ai and A2 
by uniformly randomizing their associated lengths Li, L 2 
in such a way that Li G [L® — 6 i,L^f^ + di] and L 2 G 
[^ 2 °^ — 62, +< 52 ]. In our simulations below we will use 

Aq = Ai = 1,^2 = 1.2. The transmittance T has been 
evaluated using the transfer matrix formalism Eq. Q. 

Erom the Anderson theory of localization we expect 
that for long enough tubes (and/or strong disorder) the 
transmittance will decay exponentially with the size of 
the disordered sample L (associated with N ^ L/ (L^°^ -|- 
^ 2 *^^) number of scattering units). This exponential decay 
is best described by the rescaled localization length ^ = 
+ ^ 2 °^) which is defined as: 

r(o) , ,( 0 ) 

lim ^ ^ (ln(r)). (9) 

L—^oo Jj 

where (■ ■ ■) indicates an averaging over disorder realiza¬ 
tions. All results have been averaged over more than 200 
different disorder realizations. 

The dependence of ^ on the scaled fluid velocity 
^(0) _ y{o)jf,^ fQj. various disordered strengths 61,62, 
mean lengths L^i\ and different frequencies is re¬ 
ported in the main part of Fig. Our results indicate 


where the subscript {N — 1) designate the transmis¬ 
sion and reflection amplitudes for the {N — 1) chain, 
while superscrips (1) (and Li) refer to the added scat¬ 


tering unit. Writing r^^^ = \/“^e and r 

{N-l) 


7 


(TV-l) 

R 


'yR'p ’■« ', yields the following expression for 

the logarithm of the total transmittance 




InTi 


(AT) 


= InT; 


( 1 ) 


InTi 


(AT-l) 


( 11 ) 


- In 


1 -b cos$ 


where $ = (k( + -|- cp^a) + (p^(iv-i) and is 

given by Eq. Q with the substitution {Ad,Ld,(pd) —t 
{A 2 , L 2 ,<p 2 )- Eurthermore, if one assumes that random¬ 
ness in Li, L 2 is such that the associated phases (pi ^2 are 
completely randomized (i.e. uniformly distributed be¬ 
tween —TT and tt), then averaging over these phases yields 
the rescaled (inverse) localization length 


= In 


-{l + 2 ri+ ^/l + Ar]) 


( 12 ) 


Thus, in the RPA the (rescaled) localization length ^ 
(see Eq. @) depends only on the ratio A 1 /A 2 (recall 
the definition of rj, Eq. ([^). The RPA Eq. ([T2) is also 


plotted in Fig. and matches our numerical results in 
the large domain. 

The sensitivity of ^ to for small fluid velocities 
(see Fig. can be understood by considering the aver¬ 
aging over the random lengths Li, L 2 more carefully. We 
first assume that p 2 is not random (all lengths L 2 = A® 
are identical) while the randomness in Li is such that 
6 (pi = 2{u}/c)6i/{l — Pi) > TT. In this case, the propa¬ 
gation phases at the Ai sections are completely random¬ 
ized so that one can expect (InT^-^^)! « NlnT 2 where 
T 2 is given by Eq. ([^ with the obvious substitution 
{Ld,Ad,pd) —t {L 2 ,A 2 ,(p 2 ) ((•••)i indicates an average 
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40) _L r(0)^ 

versus the rescaled lead fluid velocity /co. Vari¬ 

ous lines correspond to different disordered strengths <5, mean 
corrugation width and different frequencies co. The 

horizontal black dashed line is the result of RPA Eq. (|12|). 


FIG. 2: Dimensionless localization length ^ 



FIG. 3: Variance of logarithm of transmittance versns 
its mean value (luT) for various disordered strengths (5i,52, 
mean widths L':^\ frequencies, and velocities in the 
interval [0,0.8]. All data collapse to one universal curve char¬ 
acterizing also the = 0 case (black solid line). 


over ()iii-phases). In the case of ry <C 1, the rescaled local¬ 
ization length reads 


5 « 2?7 [1 - cos((()2)] 


(13) 


which indicates that when [32 is changing, exhibits 
(aperiodic) oscillations. The m—th oscillations is com- 

pleted when /32 takes the values P 2 “ ( ^ coT^m ) 

Next we introduce randomness also into the phase 4>2, 
i.e. we assume that (j )2 = <(' 2 °^ + ^(l >2 where = 


2(a;/co)T2°V(l “ Pi) and A(l)2 = 2(a;/co)AL2/(l - PD¬ 
As long as A(j >2 ^ 1, the disorder in L 2 has only a mi¬ 
nor effect and ^ is approximately given by Eq. ( |l3| ) , with 
(j )2 replaced by its average value ■ There is though a 
small correction, related to the variance of L 2 , so that 


r' = 277 


1 - cos(4°^) 


cos((/)^°^) 


D-Pl 


■ (14) 


We should note that the above relation is not applica¬ 
ble when /I 2 —>■ 1 since in this case A (/)2 is not small. 
Moreover Eq. (141 indicates that when <C 1 then 


^ ^ « 477(a;L2 7co)^/(l—/3|)^ i.e. the localization length 
decreases as P 2 = 0/3^°^ increases. In the opposite case 
of 1 a decrease or growth of f as /32 increases 

from zero (and while P 2 ^ P) depends on the sign of 
sin(2a;L2°Vco) (see Fig. [^. This condition can be ob¬ 
tained from the expansion of cos{(jp 2 ^) for small (32 and af¬ 
ter neglecting the small contributions from the last term 


in Eq. (14|. Finally when A02 1 we recover the RPA 


also for the segments L 2 . In this case {cos{(j) 2))2 = 0 in 


Eq. (13) and we obtain the results of Eq. (121 for ry «C 1. 


We have also analyzed the fluctuations of the trans¬ 
mission for various disordered strengths , (52, system 
sizes L and fluid velocities Our numerical results 

for (Tj^j^ = ((InT)^) — (InT)^ versus (InT) are reported 
in Figj^for various values of S [0, 0.8]. We find that, 
despite the sensitivity of p to the fluid velocity (for small 
values of the variance is a universal function 

of (InT) = L/^, independent of the velocity of the fluid. 
In fact, the analysis indicates that this universal function 
is identical to the one associated with the = 0 case 
(black bold line in Fig. [^. Our results indicate strong 
fluctuations in the localized regime due to interference 
among multiple scattered sound waves. 

Conclusions -In conclusion, we study sound waves 
propagating “on top” of a stationary fluid flow, in the 
presence of disorder. It turns out that the stationary 
flow can have a significant effect on the interference pat¬ 
tern of the waves, as compared to the case when the fluid 
is at rest. This happens because the speed of sound (in 
the laboratory frame), and therefore the phases accumu¬ 
lated between scattering events, depend on the propaga¬ 
tion direction of the wave, i.e. along the flow or oppo¬ 
site to it. We find that in a broad range of parameters 
the localization length ^ is very sensitive to the flow ve¬ 
locity v^. When increases, ^ can grow, diminish or 
oscillate- depending on the precise value of other param¬ 
eters. However for large (comparable to the speed 
of sound) phases become completely randomized and ^ 
saturates at some universal value, in agreement with the 
“random phase approximation”. 
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